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Abstract: Generically, non-linear bimetric theories of gravity suffer from the same 
Boulware-Deser ghost instabihty as non-hnear theories of massive gravity. However, re- 
cently proposed theories of massive gravity have been shown to be ghost-free. These the- 
ories are formulated with respect to a flat, non-dynamical reference metric. In this work 
we show that it is possible to give dynamics to the reference metric in such a way that 
the consistency of the theory is maintained. The result is a non-linear bimetric theory of 
a massless spin-2 field interacting with a massive spin-2 field that is free of the Boulware- 
Deser ghost. To our knowledge, this is the first construction of such a ghost-free bimetric 
theory. 
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1. Introduction and summary 

Generically, non-linear bimetric theories of gravity ||^] suffer from the same Boulware-Deser 
ghost instability ||2| as non-linear theories of massive gravity. In this paper we construct 
non-linear bimetric actions that are free of the Boulware-Deser ghost. This is possible due 
to the recent progress made in constructing ghost free theories of massive gravity. 

Until recently, consistent non-linear extensions of Fierz-Pauli massive gravity |^ 
remained elusive. In a two-parameter family of massive gravity actions was proposed 
and shown to be free of the Boulware-Deser ^ ghost in the so-called decoupling limit and 
more generally at low order in perturbation theory. The absence of the Boulware-Deser 
ghost was demonstrated at the full non-linear level in [R], based on the reformulation of 



the theory given in Q. (For complementary analyses see Q, [IC]. For a review of massive 
gravity and the associated literature, see |11|.) 



To construct non-linear generalizations of the Fierz-Pauli mass term, an additional 
"reference" metric /^jy is invariably required. In the works mentioned above, this reference 
metric is taken to be flat and non-dynamical. There are many reasons one would like to 
go beyond a flat, non-dynamical /^j^. From a theoretical standpoint, promoting f^i, to a 
full dynamical metric is desirable as leads to a background-independent theory which is 
invariant under general coordinate transformations, without the introduction of Stiickelberg 
fields. 

Such theories were introduced in |||] in order to describe the interaction of gravity 
with a massive spin-2 meson. More recently, there has been renewed interest in bimetric 
theories of gravity due to their accelerating cosmological solutions (see, e.g., [|l2|). However, 
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consistent theories of interacting spin-2 fields have proved difficult to obtain |13, 14] for 
the following reason. 

Generically, bimetric theories give rise to one massive and one massless spin-2 field 
|l|, ^, 15, |l^]. In fact, it was shown in |17] that there are no consistent theories of interacting 



massless spin-2 fields. Thus theories with multiple interacting spin-2 fields will necessarily 
include a massive spin-2 field. As a result, bimetric theories are plagued with the same 



ghost problem that generically afflicts theories of massive gravity [||, 15, |l^. Thus, ghost- 
free massive gravity is a promising starting point from which to develop a consistent theory 
of interacting spin-2 fields. 

In order to consistently promote the reference metric f^^u of a massive gravity theory to 
a dynamical variable, it is first necessary to show that the massive gravity theory remains 
ghost free for a generic non-flat, non-dynamical f^u- This would guarantee that fluctuations 
of f^u will not disrupt the consistency of the sector of the theory. Such a proof was 
carried out in jl^. Next, one must introduce a kinetic term for /^j, in such a way that the 
ffiu sector of the theory is also ghost-free. This is the aim of the present work. 

We start by reviewing the recently proposed massive gravity actions formulated with 
respect to a non-dynamical reference metric. We then show that the ghost-free mass term 
for the original metric g^u can be reformulated as a ghost-free mass term for f^y. This 
implies that the correct, consistent kinetic term for f^^, is in fact the Einstein-Hilbert term. 
We next consider the minimal bimetric model. When expanded around a fixed background, 
we show that one linear combination of fluctuations of the metrics describes a massive spin- 
2 field while the other combination describes a massless spin-2 field. We go on to discuss 
the Hamiltonian and momentum constraints of the general bimetric theory. First we give 
general arguments that the full non-linear theory also describes one massive spin-2 field 
interacting with a massless spin-2 field. We then consider the minimal model in the ADM 
formulation at the full non-linear level. We solve the constraints explicitly and demonstrate 
that our general counting arguments go through. We then extend these results to the most 
general bimetric model. We end with a brief discussion of the coupling of the bimetric 
theory to matter. 

2. Non-linear massive gravity 

The most general ghost-free non-linear massive gravity theories are given by a two-parameter 
family of actions. For a fiat reference metric, these were proposed in ||5|, ^ and shown to be 
ghost free at the complete non-linear level in . We are interested in the generalization of 
these models to arbitrary /^jy which was also shown to be ghost free in at the complete 
non-linear level. For our purposes it is convenient to work with the reformulation of the 
massive action given in Q for a general ff^^,, 

4 



L n=0 



(2.1) 



where the square root of the matrix is defined such that y/g\fy/g^f = g^'^fxu- The 
Gk{\/ 9~^f) elementary symmetric polynomials of the eigenvalues A„ of the matrix 
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\J g ^f. For a 4 X 4 matrix they can be written as, 

eo{^/9^) = 1, (2.2) 
eiiVg-^f) = Ai + A2 + A3 + A4 , 



e2iV 9 V) = A1A2 + A1A3 + A1A4 + A2A3 + A2A4 + A3A4, 



esiV 9 ^f) = A1A2A3 + A1A2A4 + A1A3A4 + A2A3A4 , 
64(^5-1/) = A1A2A3A4 = det y/g-^f . 

The highest order term in the potential in (^.l]), 



V-det5/34e4(x/5^) =/34 V-det/, (2.3) 

is independent of g^i, and so does not contribute to the g/^u equations of motion. Thus when 
the reference metric f^u is taken to be non-dynamical, this term can be neglected. The 
four remaining /3„ describe four combinations of the mass of the graviton, the cosmological 
constant and two free parameters. 

In what follows, we wish to consider the case that f^u is dynamical and thus the 
potential term in (2.1) acts as a potential for f^i, as well as g^i,. Thus the term 



should be retained in the action. The resulting theory has five free parameters, including 
the graviton mass and the cosmological constants for both g^i, and f^u- 

In ref. ||l8[| , the action ( |2.1[ ) was analyzed in the ADM formulation for a general non- 
flat but non-dynamical ffj^^. In the ADM formulation, there are six potentially dynamical 



modes described by gij and their conjugate momenta tt^K It was shown in [18j and [19 
that, for these theories ( |2.lD , there is a Hamiltonian constraint along with an associated 
secondary constraint. These constraints remove one propagating mode from the theory, 
leaving only the five modes consistent with a massive spin-2 particle. Thus, for a non- 
dynamical ff^i^, the g^i, sector of theory is free of the pathological sixth mode known as the 
Boulware-Deser ghost. The existence of the constraint is due to the specific form of the 



mass term in (2.1) 



3. Bimetric gravity 

Let us now consider the possibility of adding a kinetic term for f^^. For a dynamical f^u, 
the mass term in ( p.lj ) acts as a potential for /^jy as well as for 5^,^. We want to add 
a kinetic term for that is consistent with this potential, to get a ghost-free theory. 
The given mass term does not appear to have the same form for as it does for g^j,^. 
This is not automatically calamitous, however, as we are not constrained to write only the 
Einstein-Hilbert kinetic term for /^j^. In general, kinetic terms for /^^^ can be constructed 
using both metrics g^^ and /^^. It could be that there exist multiple kinetic terms that 
give rise to a ghost-free theory in the f^^ sector. We wish to find one such kinetic term. 

To determine what valid kinetic term we should introduce for f^^, it is first helpful to 
rewrite the above mass term so that it more closely resembles a potential for /^j^. Using 
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t/— det g = \J — det / \l det / ^(7, we have, 

4 4 

y^dit^5^/3nen(\/5^) = y^dit7e4(A/rt) ^/3„e„(\/^). (3.1) 



n=0 n=0 



We would also like to express the en{_\J g~^f) in terms of the inverse matrix f~^g- Note 
that the elementary symmetric polynomials of f~^g can be easily expressed in terms of 
the eigenvalues of g~^f, 



eo 



iVF^g) = 1, (3.2) 



/ ^r^^ 1111 
Ai A2 A3 A4 

e2(\/7^) - ^ + ^ + + + + ^ 
A1A2 A1A3 A1A4 A2A3 A2A4 A3A4 

1111 



e4(^/rt) 



A1A2A3 A1A2A4 A1A3A4 A2A3A4 
1 



A1A2A3A4 

Comparing ( |2.2| ) and ( |3.2D , it is evident that the en{\/ g^^ f) can be expressed entirely in 
terms of the en{\fT^), 



e.(v^) = ''-fiJ^'^ ■ (3.3) 



Then it follows that the mass term ( |3.1| ) can be re-expressed as, 

4 4 

\/^d^^/3nen(\/^) = y^d^^/5ne4-„(yrt). (3.4) 

Now it is apparent that the mass term does in fact have the same ghost-free structure for 
f^u as for g^y, only with different coefficients. This implies that a valid kinetic term for 
ffj_y is the standard Einstein-Hilbert kinetic term after all. Introducing such dynamics for 
f^y will give a theory that we expect to be free of the Boulware-Deser ghost in both the 
g^u and f^^ sector, 

S = j d^xyj- det g i?'"' + j d^x^/- det / i?*-^' 

+2rr?Mls / (fx^ - det gY, Pn eniVg^f) ■ (3.5) 

n=0 

Here -R*^' denotes the scalar curvature for g^i, and ii*^' denotes the scalar curvature for /^j,. 
We have introduced different Planck masses for the two metrics and defined an "effective" 
Planck mass as. 
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Given the arguments of ||5|, |6|, 0, |T^, the mass term ( |3.1| ) is the unique interaction term 
that avoids the Boulware-Deser ghost in the (^-sector with a non-dynamical /. Here we 
have shown that it also avoids the ghost in the /-sector with a non-dynamical g. 

This structure suggests that the bimetric action is free of the pathological Boulware- 
Deser modes. However, when both g and / are dynamical, it is not yet obvious how the 
mass term avoids the two Boulware-Deser ghosts simultaneously. In addition, the above 
arguments are not sufficient to determine the total number of propagating modes of the 
theory at the non-linear level. We examine these issues in the following sections. 

As for the uniqueness of the kinetic term for f^y, it remains to be seen if there exist 
other such terms, or if consistent kinetic mixings can be introduced between g^u and f^^. 
However, given the theorem of [0], such terms must vanish in the limit m — )• 0, while higher 
derivative mixings can be neglected compared to the potential in the long-wavelength limit 
ll^. Consequently, we consider only the action given by ( |3.5D in what follows. 



4. Spectrum of the linear theory 



Bimetric theories invariably describe one massive and one massless spin-2 field [||, T^, 16 1. 
This is apparent for any theory which reduces to the Fierz-Pauli theory ^ at the linear 
level. There, the fluctuation g^j^^ — f^y acquires a mass as this is the only combination of 
metrics that appears in the potential. 

To see this explicitly, let us consider small perturbations in the minimal massive model 
introduced in [^. The interaction term of the minimal model is given by 

2m2M2j j d^x^J-deig (s - tr 7^ + det , (4.1) 

corresponding to the coefficients, 

/3o = 3, /3i = -l, /32 = 0, /33 = 0, /34 = 1. (4.2) 



Let us expand both g^y and f^y around the same fixed background g^y, 



To second order in perturbations, the minimal model gives 

S = j (fx ( V'^'^'"''^"/? + l^^ui^'^'^lap) 



(4.3) 



(4.4) 



d^x 



Mr, M 



Mr, M 



Here £'^'^"/^ denotes the usual Einstein-Hilbert kinetic operator and indices are raised and 
lowered with respect to the background metric. 

To diagonalize, we perform the change of variables 



1 



1 



1 



1 



Mf + Mg ' 

" Mg 



Mf ' 



(4.5) 
(4.6) 
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The action (0) becomes 

S = Jd'x [u^J^^^^u^p + v^.J^'^'^^v^p) -^Jd'x [v^^v^y - v^^v\) . (4.7) 

Thus at the hnearized level, we see that the minimal action describes one massless spin-2 
particle Ufj^iy and one massive spin-2 particle Vf^^^ with mass m. 
Note that, 

for Mf > Mg : u^^ 1^^ and v^u h^y , 
for Mg > Mf : u^y h^u and v^u -l/iu ■ 

When one Planck mass greatly exceeds the other, it is always the massless particle that has 
the larger Planck mass. Thus it is possible to send one Planck mass to infinity so that the 
massless particle decouples from the theory. One is left with a formally covariant theory 
with a single massive spin-2 field propagating in a fixed background. 



5. Constraint analysis: minimal theory 

Beyond the linear level, one can determine the number of degrees of freedom of the bimetric 
theory by doing an analysis of the Hamiltonian and momentum constraints in the ADM 



formulation [20]. Such an analysis will determine whether or not the Boulware-Deser ghosts 
avoided at the linear level reappear in the non-linear theory. Let and Ni denote the lapse 
and shift functions of the metric g^u while L and Lj denote the lapse and shift functions 
of the metric Accordingly we define, 

= (-5°°)"'/' , Ni = goi, l^j=9ij, (5.1) 
L = {-f°r'^\ L, = U, %j=fij. (5.2) 

The canonically conjugate momenta for jij and ^fij are given by vr*-^ and p*-^ respectively. 
In the bimetric theory ( |3.5| ) the lapse and shift variables appear without time derivatives 
and are thus non-dynamical. The remaining six gij, six fij and their conjugate momenta 
constitute 24 phase-space degrees of freedom, i.e., 12 potentially propagating modes, a 
propagating mode referring to a pair of conjugate variables. 

Before before proceeding with the constraint analysis of the minimal model we note 
that, based on the arguments of the previous sections, we expect the counting of degrees 
of freedom to be appropriate for a massless spin-2 field (two modes) plus a massive spin-2 
field (five modes). Thus, we need ten total constraints and gauge invariances to leave us 
with 14 degrees of freedom, corresponding to seven total propagating modes. The general 
covariance should be responsible for removing eight degrees of freedom. Four degrees of 
freedom will be removed by gauge fixing while the other four will be removed via the 
associated Hamiltonian and momentum constraints. Indeed, we will show below that for 
the minimal theory, one set of lapse and shift variables appear linearly, thus their equations 
of motion act as constraints on the remaining variables. In addition, we expect there to be a 
single constraint and associated secondary constraint, as was found in |18| and |19]. These 
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remove an additional propagating mode. Below we will also demonstrate the existence of 
this constraint in the minimal theory. Thus, out of 24 initial degrees of freedom we are left 
with 14, consistent with one massive spin-2 field and one massless spin-2 field. 

Let us again be more explicit by studying the constraints of the minimal model ( [4.1] ) . 
This action is highly non-linear in the lapse and shift variables. However, following ||l^ we 
introduce three new variables n* which are functions of the lapse and shift variables and 
spatial metrics, 

-V = {L5\ + ND\)n^ . (5.3) 

The matrix D = D\ is given by, 



(5.4) 



where 



fmj) ■ 



(5.5) 



Significantly, D is independent of the lapses N and L, and the shifts A^* and L*. Let us 
also define 



X 



1 



,k a 



n 



Jkm)n 



(5.6) 



We can use the above variables to replace the shift N"^ in favor of the n*. The minimal 
action then becomes linear in the lapses N and L, and the shift L*, 



+L 
+N 



M^R^'^'+M^Din'^R!^' 



2m^M^s{^/d^V^ - Vdet 'f) 



2m^M^s{^/det^^/xD\ - S^/detj)] . (5.7) 



What's more, the equations of motion for the can be shown to be independent of N 
1 18]. Thus they can be used to fix in terms of the remaining degrees of freedom. 



n 



-M^Cf-'Y^R'; 



is) 
'3 



Am^M^Q det 7 + M^R['^ (i-^ f^R\' 



The equation of motion is 



r2 d's^ 



xD\ 



(5.8) 



(5.9) 



On substituting for n* we can use this equation to constrain the vr*-^ and 'jij. Along with 
the associated secondary constraint found in |1£], this eliminates one propagating mode. 

A crucial point is that the ghost mode of 7ij that is eliminated through the Hamiltonian 
constraint ( |5.9D remains independent of the lapse L. Only then does the term Tr'^^dt'jij in 
the Lagrangian (^3)' when evaluated on the constraint surface, remain independent of L. 
To see this, note that n*, ^/x and D are independent of the lapse L. Thus, the Hamiltonian 
constraint ( ^.91) and the associated secondary constraint are independent of L as well. Using 
the constraint ( |5.9| ), the would-be ghost mode of the (7-sector can be expressed in terms of 
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the remaining independent variables. It follows that the expression for the would-be ghost 
mode does not involve L either. 

Now we use equations (|5.80 and ( |5.9D to eliminate the lapse and the new shift-like 
variables n* of the (/-metric from the Lagrangian (5.7). The resulting expression remains 
linear in both the lapse L and shift L* of the /-metric, 



(5.10) 



' det 7 



Thus U and L are Lagrange multipliers. Their equations of motion act as four constraints 
on the six canonical pairs {p^-^ ,^fij) and the remaining five independent pairs of (n^^jjij), 



R^i'^+ M] R^f=0, 



(5.11) 
(5.12) 



Note that using ( |5.1l| ), we can rewrite ( ^.12 ) in terms of the /^jy variables i?^^^ as 



Mj 



2m2M2j 



o(/) 



' det 7 



M^ 



Am^M^^ 



(5.13) 



The situation now resembles that of general relativity. These four constraints along 
with four coordinate conditions can be used to eliminate another four pairs of canonical 
variables, including the ghost in the /-sector. Thus we are indeed left with seven propagat- 
ing modes, consistent with one massive and one massless spin-2 field. Roughly speaking, 



(|5.12D can be used to eliminate the Boulware-Deser mode of the /-sector while the con- 
straints ( 5.11| ) eliminate one combination of the helicity-1 and helicity-0 parts of g^y and 
f^u- The surviving combination, together with the corresponding combination of helicity-2 
modes of g^y and /^i^, will make up the massive state. After fixing the coordinate condi- 
tions, the L and U are determined by four of the ^fij and p^^ equations of motion. Another 
four of these equations should be trivially satisfied as a consequence of Bianchi identities. 

Note that we have introduced our new shift-like variables ( |5.3D in such a way that an 
asymmetry arises between the ^-sector and the /-sector in the ADM analysis. Namely, 
there appears to be only one constraint arising from the (/-sector of the theory while there 
are four constraints in the /-sector. Let us emphasize that this apparent asymmetry is 
purely a result of the choice of variables (|5.3D and that no fundamental asymmetry exists 



between these two sectors, as is evident from the constraint equations ( ^.91) , ( 5.11 ) and 



0. 



6. Constraint analysis: full theory 



Let us now extend our analysis to the full bimetric theory (3^). We write this action in 



terms of the ADM variables and again replace the shift N"^ with n* as defined in (5.3). The 
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full Lagrangian is then linear in the lapses N and L, and the shift U, 

c = Ml Tr'^daij + M] p'^dt%j -no + NC, 

where T-Lq and C stand for 



C = mIr^'"+ M^Ri'Din^ + 2m^M^ff^/det^V . 
U and V in the above expressions are defined as 



n 



+133 {D\ %jD^,n>' - D\ %,D\n') + \^\d\D' ^ - D^^D\] 



+/34Vdet^/\/dit^, 



Varying the Lagrangian with respect to N gives the equation of motion, 

C = 0, 

while the equation of motion is given by 



0. 



(6.1) 

(6.2) 
(6.3) 

(6.4) 
(6.5) 



(6.6) 



(6.7) 



This equation is independent of the lapses and L, and the shift L* and can in principle 
be used to determine the in terms of the dynamical variables jij and vr*-^. Thus the A'^ 
equation (^^) becomes a constraint on the dynamical variables, reducing the 12 degrees of 
freedom to 11. This is the Hamiltonian constraint. 

In addition, because n*, -y/x and D are independent of the lapse L and shift U, the 
Lagrangian ( |6.lD remains linear in L and U after the n* have been eliminated in favor of 
7jj and tt^^ . Thus the lapse and shift remain Lagrange multipliers, enforcing the correct 
number of constraints, along with the general coordinate invariances, and primary and 
secondary constraints, for one massive and one massless spin-2 field. 



7. Coupling to matter 

The coupling of matter to the metrics 5^,^ and f^i, must not violate the constraints found 
above. Hence the coupling should be linear in the lapse and shift functions of both g^i, and 
ff^^. The minimal coupling of General Relativity to matter, ^/— det g Cm, automatically 
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satisfies this requirement. If this were not the case, the constraints of even massless GR 
would be violated. Thus the simplest allowed matter couplings are of the form, 

V^d^jOiig^,,^^), V-det//:2(^.,0. (7.1) 

Here the ^"^'^ represent any matter field. 

If both g^i, and /^j^ are allowed to couple to the same matter fields, the particle 
geodesic equations will be modified, generically violating the equivalence principle. Thus 
in many cases such couplings could be easily ruled out experimentally. This problem could 
be avoided if matter couples to an effective metric gfj,i, constructed out of and /^jy, as 

^/-detg C3{gf,,,(l>^). (7.2) 

However, such a coupling will generically violate the constraints. Thus the allowed cou- 
plings to matter in the bimetric theory are severely restricted. The observational implica- 
tions of these couplings and possible classical solutions will be discussed in future work. 
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